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1.  INTRODUCTION 

A  team  of  ten  Ph.D.  scientists  and  several  graduate  students  has  been  assembled  at 
USU  in  order  to  establish  a  ‘Center  of  Excellence  in  Theory  and  Analysis  of  the  Geo- 
Plasma  Environment’  The  USU  scientists  work  in  close  collaboration  with  colleagues  at 
the  Phillips  Laboratory  (PL)  in  Bedford,  Massachusetts,  on  a  number  of  problems  that  are 
relevant  to  Air  Force  systems,  including  OTH  radars,  communications,  and  orbiting  space 
structures.  The  overall  goal  of  the  research  is  to  obtain  a  better  understanding  of  the  basic 
chemical  and  physical  processes  operating  in  the  geoplasma  environment,  including  the 
ionosphere,  thermosphere,  and  magnetosphere.  Some  of  the  more  specific  goals  are  as 
follows: 

1 .  Study  the  production,  transport  and  decay  of  ionospheric  density  structures. 
This  will  include  studies  of  Ae  main  electron  density  trough,  sun-aligned  arcs, 
and  plasma  blobs. 

2 .  Assist  in  the  development  of  better  operational  ionospheric  models  for  the  Air 
Force.  This  includes  running  our  numerical  model  of  the  global  ionosphere  for 
numerous  geophysical  cases  so  that  the  results  can  be  incorporated  into  the 
AWS  Ionospheric  Specification  Model. 

3 .  Conduct  model/data  con:q)arisons  in  order  to  validate  the  ionosphoic  models. 

4.  Develop  a  ray  tracing  code  for  our  global  ionospheric  model  so  that  we  can 
support  PL’s  OTH  raito  campaigns. 

5 .  Explore  methods  of  incorporating  the  PL  particle  precipitation  energy  index  in 
bodi  magnetospheric  convection  and  particle  precipitation  models,  which  are 
inputs  to  the  ionospheric  and  thermospheric  specification  models. 

6.  Study  plasma  convection  characteristics  at  high  latitudes  with  emphasis  on 
northward  interplanetary  magnetic  field  (IMF)  conditions. 

7 .  Study  magnetosphere-ionosplKre  coupling  problems,  including  wave  excitation 
associated  vtith  the  polar  wind  and  energetic  ion  ou^ows.  We  will  also  study 
the  electrical  coupling  and  currents  in  the  magnetosphere-ionosphere  system. 

8.  Continue  the  development  of  a  high-resolution  multi-species  thermospheric 
circulation  model.  'Dus  is  relevant  to  AWS’s  interest  in  developing  a  real-time 
Thermospheric  Specification  Model. 

9.  Study  plasmasphere  refilling  characteristics  associated  with  magnetic  storms 
and  substorms. 

10.  Study  certain  spacecraft-environment  interaction  problems,  particularly  those 
relat^  to  solar  cell  changing,  high-voltage  power  sources,  artificial  plasma 
cloud  expansion,  and  outgassing  from  large  structures  orbiting  at  ionospheric 
altitudes. 
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2.  USU-PmLUPS  LABORATORY  COLLABORATIVE  EFFORTS 

At  the  present  time,  we  have  ongoing  collaborative  efforts  with  several  scientists  at 
the  Phillips  Laboratory  in  Bedford,  Massachusetts.  In  the  paragraphs  that  follow,  we 
briefly  describe  the  progress  made  during  die  last  year. 

2.1.  D.  Anderson  (LIS) 

We  helped  D.  Anderson  develop  a  High  Latitude  Ionospheric  Specification  Model 
for  Air  Weather  Services.  Our  original  effort  involved  running  our  numerical  ionospheric 
model  for  both  the  northern  and  southern  hemispheres  for  a  range  of  solar  cycle,  seasonal, 
magnetic  activity,  and  IMF  conditions.  A  total  of  108  ionospheric  simulations  were  run. 
R.  Daniell  of  (Computational  Physics  Incorporated  ((CPI)  *fit’  our  calculated  densi^ 
distributions  as  a  function  of  altitude,  latitucte,  longitude,  season,  solar  cycle,  magnetic 
activity,  and  univorsal  time  using  relatively  sitiqile  functions.  These  background  density 
distributions  were  then  incorporated  into  his  Parameterized  Real-Time  Ionospheric 
Specification  Model  (PRISM).  During  the  last  year,  D.  Anderson  requested  that  we 
determine  the  accuracy  of  die  fitting  procedure  and  D.  (Crain  in  our  group  undotook  this 
effort  D.  (Crain  found  that  the  paramererized  model  developed  by  R.  Daniell  reproduced 
the  original  density  distributions  very  well  and  that  the  parameterized  version  of  our 
numerical  ionosphoic  model  was  much  more  rejnesentative  of  the  high-latitude  ionosidiere 
than  existing  m^els. 

At  the  request  of  D.  Anderson,  J.  Sojka  is  spearheading  an  inter-model  comparison 
of  the  noon  sector  mid-latitude  F  region.  This  inter-model  comparison,  PRIMO,  is  a 
CEDAR  activity  but  is,  to  a  large  extent,  the  outgrowth  of  USU-PL  interactions  in 
ionospheric  model  and  data  comparisons.  The  problem  being  focussed  on  in  this  research 
area  is  the  difficulty  during  solar  maximum  coriditions  of  reproducing  the  observed  NnjFi 
values.  The  models  are  generally  too  low.  A  presentation  of  these  results  was  made  at  the 
(CEDAR  1991  meeting. 

2.2.  S.  Basu  (US) 

In  collaboration  with  S.  Basu  and  C.  Valladares  (Boston  College)  we  studied  the 
properties  of  sun-aligned  arcs.  In  one  effort,  the  emphasis  was  on  the  ionospheric 
modification  caused  by  a  typical  sun-aligned  polar  cap  arc.  The  radar,  optical,  and  satellite 
data  pertaining  to  these  arcs  that  were  available  at  Phillips  Laboratory  were  used  to  define 
the  arc  features  and  our  ionospheric  model  was  then  run  using  these  parameters  to  elucidate 
the  ionospheric  response  to  the  arc.  A  paper  describing  this  work  is  currently  being 
written. 

2.3.  H.  Carlson  (LIS)  -> 

At  the  beginning  of  this  year’s  effort,  H.  Carlson  emphasized  the  need  to 
understand  the  electrodynamics  of  a  sun-aligned  polar  cap  arc.  He  was  interested  in 
understanding  how  the  system  self-consistently  evolves  from  a  fairly  uniform  ionosphere 
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to  a  polar  cap  arc  situation.  He  was  also  interested  in  determining  what  contribution  the  F 
region  ionization  makes  to  the  overall  conductivity  of  the  sun-aligned  arc.  In  response  to 
these  interests,  L.  Zhu  initiated  the  development  of  a  time-dependent  nK)del  of  a  sun- 
aligned  arc  in  which  the  electrodynamics  are  treated  self-consistently  in  the  frame  of  the 
magnetosphere-ionosphere  system.  The  model  development  was  conq>leted  during  the  last 
year  and  die  model  is  currmtly  being  used  to  study  the  tenoral  development  and  evolution 
of  an  arc  driven  magnetospheric  processes. 

2.4.  W'.  Demg  (PHG) -4 

The  very  strong  geomagnetic  storm  of  13  March  1989  has  been  tiie  focus  of  a  close 
collaboration  with  Dr.  W.  Doug.  The  first  drafi  of  the  publication  is  well  in  hand.  The 
extensive  DMSP  data  base  for  this  March  1989  period  has  been  interpreted  by  W.  Denig 
and  used  both  as  an  input  data  set  to  our  global  iono^heric  nxxiel,  and  as  an  in  situ  plasma 
density  reference  data  set  to  compare  against  the  model  results.  The  pre-storm  observation 
and  model  agreement  is  good.  During  the  storm  phase  the  comparisons  are  very  good. 
Both  DMSP  and  TEC  data  were  used  for  die  comparisons. 

2.5.  J.KlobucharQJS)-^ 

In  conjunction  with  J.  Klobuchar  at  Phillips  Laboratory  and  P.  Doherty  at  Boston 
College,  we  began  a  study  of  the  TEC  variation  observed  at  various  high-latitude  stations 
conpared  to  diose  predict  by  our  ionospheric  model.  Specifically,  TEC  values  predidted 
by  our  time-dependent  ionospheric  model  were  compart  with  the  diurnal  variations  of 
TiSC  measured  at  the  Goose  Bay,  Hamilton  Bay,  and  Poker  Flat  observing  sites  and,  in 
general,  very  good  agreement  between  model  and  data  was  achieved.  A  paper  on  this 
con^arison  is  currendy  being  written. 

2.6.  Drs.  D.  Hardy  and  S.  Gussenhoven  (PHP)  -> 

Captain  T.  Frocminckx  con^leted  his  Masters  Degree  at  USU  in  the  spring  of  1991 
[Frooninckx,  1991].  His  thesis  concerned  the  spacecraft  charging  experience  on  the 
DMSP  satellites.  Captain  I^ooninckx  worked  with  scientists  at  the  Phillips  Laboratory, 
specifically,  Drs.  D.  Hardy  and  S.  Gussenhoven.  As  a  result  of  this  work,  a  strong  solar 
cycle  dependence  in  the  spacecraft  charging  was  postulated  and  verified  using  the  DMSP 
data.  Two  papers  were  written  and  submitted,  one  has  now  been  published  (papers  30  & 
35). 

2.7.  J.  Buchau  (LIS)  -> 

In  support  of  J.  Birchau’s  OTH  radar  work.  Captain  P.  Citrone  has  interfaced  our 
time-dependent,  high-resolution  ionospheric  model  with  the  Jones  and  Stephenson  [1975] 
ray-tracing  program  to  examine  the  effects  of  3-D  electron  density  gradients  on  ray 
propagation.  We  have  carried  out  initial  validation  tests  of  the  ray  paths  through  this 
ionosphere.  These  tests  indicated  an  azimuthal  position  bias  as  the  ray  paths  are  bent  by 
h(»iz(mtal  electrcm  density  gradients.  Work  is  underway  to  simulate  iono^heric  conditions 
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for  a  PLrOTH  campaign.  The  potential  value  of  this  physical  model  is  that  it  can  generate 
physically  realistic  ionospheric  conditions,  which  can  be  used  as  a  refraence  for  ray-tradng 
simulations  on  a  specific  day,  time,  and  location.  In  situ  data  can  then  be  used  to  modify  or 
add  fine  structure  as  appropriate.  Or,  in  the  absence  of  data,  parametric  sensitivity  studies 
can  readily  be  carried  out  For  his  M.S.  thesis,  P.  Gtrone  studied  the  effect  of  the  trough- 
auroral  boundary  on  the  propagation  and  bending  of  HF  radio  waves. 

2.8.  D.  Cooke  (PHK)  -> 

We  are  studying  a  range  of  spacecraft-environment  interaction  problems  that  are 
relevant  to  D.  Cooke’s  interest  These  inclucte  outgassing  from  the  Space  Shuttle,  the 
expansion  of  artificial  plasma  clouds,  solar  cell  charging  processes,  and  high-voltage 
sphoes  in  the  iono^here.  Several  papers  have  been  written  on  these  topics  during  the  last 
year  and  these  are  Mefly  described  in  the  next  section  and  listed  in  the  URI  Publications 
section. 

2.9.  y.  Whalen  (LIS) 

J.  Sojka  met  with  J.  Whalen  during  his  July  19, 1991  visit  to  Hanscom  Air  Force 
Base.  They  discussed  J.  Whalen’s  recent  work  on  die  extensive  IGY  ionosonde  data  base, 
which  covers  both  the  northern  and  southern  hemispheres.  They  discussed,  at  some 
length,  his  equinox  data,  which  show  the  same  Kp  dependence  for  the  afternoon  sector 
troughs  in  the  northern  and  southern  hemispheres.  Indeed,  the  conjugacy  is  remaricably 
strong  considering  the  various  shoitcomings  in  global  coverage  and  available  days  in  the 
equinox  periods  during  IGY.  We  discussed  a  joint  follow-up  to  this  work  using  the  USU 
nordiem  and  southern  hemisphere  conqiuter  simulations  that  were  tun  for  D.  Andoson  (the 
108  simulations). 

J.  Whalen  has  also  been  looking  at  the  stations  in  the  polar  ctqi  which,  during 
winter  conditions,  see  a  peak  densify  at  magnetic  noon  rather  dian  geographic  (solar)  noon. 
This  maximum  is  a  convection  feature.  He  has  data  from  up  to  5  stations  at  different 
Imgitudes  but  at  the  same  magnetic  latitude  (~7S*)  in  the  soutiiem  hemisphoe.  Each  has  a 
distinctly  different  solar  -  magnetic  noon  offset  and  each  shows  tiie  magnetic  noon  peaking. 
These  data  relate  to  the  phenomenon  referred  to  as  the  “tongue  of  ionization.’’  This  work  is 
at  an  early  stage,  but  we  are  interested  in  following  up  with  an  ionospheric  model  study  and 
a  subsequent  comparison  with  his  data  to  see  if  our  ionospheric  model  properly  describes 
the  ‘tongue’  feature. 
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3.  AWS  PERSONNEL 

Four  Air  Force  personnel  from  Air  Weather  Services  have  completed  their  M.S. 
Theses  during  the  last  two  years.  Their  names,  theses  titles,  major  professors,  and 
completion  years  are  listed  below. 

1.  Gary  R.  Huffmes 

Tide:  Using  the  USU  Ionospheric  Model  to  Predict  Radio  Propagation 
Through  a  Simulated  lono^here 
Date:  1990 

Major  Professor  Jan  J.  Sojka 

2.  Thomas  B.  Frooninckx 

Tide:  Higb-ladtude  Spacecraft  Charging  in  Low-Earth  Polar  Orbit 
Date:  1991 

Major  Professor.  Jan  J.  Sojka 

3.  David  R.  Payne 

Tide:  Electron  Heating  and  Ion  Production  Rates  in  Auroral  and  Sun- 
Aligned  Arcs 
Date:  1991 

Major  Professor  Robert  W.  Schunk 

4.  Peter  J.  Qtrone 

Tide:  The  Effect  of  Electron-Doisity  Gradients  on  Propagation  of  Radio 
Waves  in  die  Mid-Latitude  Trou^ 

Date:  1991 

Major  Professor  Jan  J.  Sojka 

In  addition,  two  new  Air  Force  students,  W.  Cade  and  M.  Loveless,  have  recently  joined 
our  group  and  will  begin  woridng  on  research  topics  for  dieir  M.S.  theses  in  the  Fall. 
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4.  SCffiNTMCACCOMPUSHMENTS 

In  addition  to  the  USU-PL  collaborative  efforts  that  are  in  progress,  we  have  also 
completed  a  number  of  additional  studies.  During  the  second  year  of  the  URI  program,  17 
scientific  papers  have  been  submitted  for  publication,  21  presentations  have  been  given  at 
both  national  and  international  meetings,  and  there  have  been  several  trips  to  the  Phillips 
Laboratory  in  Boston,  Massachusetts  and  to  AWS  in  order  to  coordinate  activities.  Lists  of 
the  URI  personnel,  publications,  presentations,  and  trips  ate  attached. 

In  the  subsections  that  follow,  we  briefly  highlight  some  of  the  scientific 
accomplishments  and  discuss  the  status  of  work  in  some  of  the  inq>ortant  areas.  The 
emphasis  during  the  second  year  of  our  URI  program  was  on  spacecraft-environment 
interaction  studies,  including  spacecraft  charging  (papers  30  &  35),  high  voltage  ^heres 
(papers  21, 23,  &  26),  arc  discharging  of  solar  arrays  (paper  27),  and  artificial  plasma 
cloud  expansion  (paper  31). 


4.1  Spacecrcft  Charging 

Spacecraft  charging  of  Air  Force  satellites  can  lead  to  serious  operational  anomalies 
and,  hence,  this  area  of  research  has  direct  and  practical  applications.  Most  of  tiie  charging 
occurs  along  geosynchronous  orbits  as  the  satellites  encounter  the  energetic  particles  in  the 
radiation  belts,  but  recent  experimental  evidence  has  shown  that  satellites  in  low-earth-orbit 
can  also  charge  to  significant  potentials.  In  particular,  experimental  evidence  clearly 
indicates  that  Defense  Meteorological  Satellite  Bogram  (DMSP)  polar  orbiting  spacecraft  at 
840  km  can  develop  electric  potentials  as  severe  as  -1430  V  while  at  high  magnetic 
latitudes.  To  explore  this  charging  region,  an  analysis  of  DMSP  F6,  F7,  F8,  and  F9 
satellite  precipitating  particle  and  ambient  plasma  measurements  taken  during  periods  of 
high,  medium,  and  low  solar  flux  was  performed.  One  hundred  eighty-four  charging 
events  ranging  from  -46  to  -1430  V  were  identified,  and  an  extreme  solar  cycle 
dependence  was  found,  as  charging  is  most  frequent  and  severe  during  solar  minimum. 
Satellite  measuronents  and  time-dqwndent  ionos^eric  model  (TDIM)  output  were  used  to 
determine  the  cause  of  the  solar  cycle  dependence  and  to  characterize  the  environments 
which  both  generate  and  inhibit  tiiese  potentials.  The  electron  precipitation  associated  witii 
various  DMSP  charging  levels  was  analyzed;  it  was  suggested  that  precipitating  electrons 
as  low  as  2  to  3  keV  may  contribute  to  charging,  although  higher-energy  electrons  make 
greater  contributions.  Secondary  electron  production  due  to  incident  electrons  below  1  keV 
was  shown  to  inhibit  charging.  The  energetic  electron  fluxes  shown  to  generate  charging 
did  not  vary  significantly  over  the  solar  cycle.  Instead,  DMSP  ambient  plasma  data  and 
TDIM  generated  results  indicated  that  a  variation  in  the  thermal  plasma  density  of  over  1  or 
more  ordoa  of  magnitude  was  the  cause  of  the  solar  cycle  dependence,  and  an  ambient 
plasma  density  of  less  than  10^  cmr^  was  found  necessary  for  significant  negative  charging 
^  100  V)  to  occur. 
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4.2  Artificial  Plasma  Chud  Expansion 

We  conducted  three-dimensional  simulations  of  the  expansion  of  artificial  plasma 
clouds  in  the  ionosphov  at  low-earth-orbit  (LEO)  altitudes  (paper  31).  Such  clouds  can  be 
created  by  spacecraft  outgassing  followed  by  ionization  or  by  the  direct  injection  of  a 
neutral  gas  into  the  ionosphere.  In  the  study,  a  3-dimensional,  time-dependent  fluid  model 
was  used  to  study  the  ionization  and  plasma  expansion  characteristics  of  barium  clouds. 
Neutral  gas  clouis  with  a  total  mass  of  1  kg  were  released  with  a  spherical  Gaussian 
doisity  distribution  and  the  subsequent  photoionization  1^  solar  UV  radiation  and  ion  cloud 
expansion  were  modelled.  Three  cases  were  considn^;  A  cloud  without  a  directional 
velocity;  a  clmid  widi  an  initial  velocity  of  5  km/s  across  the  B  field;  and  a  cloud  with  initial 
velocity  components  of  2  km/s  both  along  and  across  the  B  field.  For  the  cloud  without  a 
directional  velocity^  the  ionization  occurs  in  a  spherical  volume.  The  resulting  Ba^  cloud 
expands  along  the  B  field  and  the  Ba***  density  ^stribudon  gradually  becomes  ellipsoidal 
from  the  inner  to  the  outer  parts  of  the  Ba'*'  cloud.  The  electrostatic  snowplow  effect 
associated  with  the  expanding  Ba'*'  cloud  creates  a  hole  in  the  G**  background  at  the  center 
of  the  Ba'*'  cloud  and  creates  two  0+  density  bumps  on  Ae  two  sides  of  the  Ba***  cloud.  For 
the  neutral  gas  release  with  an  initial  velocity  across  the  B  fiekl,  the  resulting  Ba'*'  cloud  has 
a  cennet-like  density  distribution  at  eariy  times.  Eventually,  because  of  the  expansion  along 
the  B  field,  the  Ba***  cloud  becomes  sheet-like.  Again,  there  are  two  Q**  density  enhanced 
regions  and  one  O**  density  depletion  region  in  the  background  ionosphere.  For  these  two 
cases,  although  there  are  C)+  density  depletions,  tiiere  are  no  electron  density  (total  plasma) 
depletions.  J^en  the  neutral  gas  cloud  has  initial  velocity  components  both  along  and 
across  the  B  fiield,  the  situation  is  quite  different.  The  resulting  Ba'*’  cloud  has  a 
complicated  density  distribution.  The  Ba'*'  snowplow  effect  creates  an  G**  density  hole  on 
one  side  and  an  G**  density  bump  on  the  other  side  of  the  expanding  Ba^  cloud.  There  is  a 
net  plasma  depletion  on  the  side  opposite  to  tiie  Ba'*'  cloud  motion  along  the  B  field. 


4.3  Solar  Cell  Operation 

We  published  one  paper  that  dealt  with  the  interaction  of  solar  cell  modules  with  the 
ionospheric  plasma  (paper  27).  For  future  space  missions,  elevated  voltages  of  hundreds 
of  volts  are  necessary  to  minimize  mass  requirements  and  resistive  losses  of  orbiting  solar 
arrays.  However,  these  elevated  voltages  (both  negative  and  positive)  can  yield  hazardous 
conditions,  including  anomalous  current  surges,  arcing,  and  continuous  power  drains. 
Most  of  the  work  done  in  the  past  dealt  with  'positive*  voltages,  but  recently  there  has  been 
an  increased  interest  in  ‘negative’  solar  cells.  Therefore,  we  examined  both  experimental 
data  and  particle-in-cell  (PIQ  computer  simulaticms  in  (Hder  to  determine  in  what  region  of 
the  solar  array  is  arcing  likely  to  occur.  In  genoal,  it  could  occur  at  the  interconnector- 
dielectric  interface,  solar  cell  edge,  or  gap  regions.  However,  for  ‘negative’  voltages  we 
found,  contrary  to  what  was  found  for  positive  voltages,  that  arcing  occurs  around  the  gap 
region  between  adjacent  solar  cells  and  also  at  the  solar  cell  edge  -  substrate  interface,  but 
not  around  the  negatively  biased  interconnectors. 


4.4  High-Voltage  Spheres 
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Three  of  our  papers  were  concerned  with  the  interaction  of  a  high-voltage  sphere 
with  die  LJEO  plasma  environment  (papers  21, 23.  &  26).  This  woik  was  motivated  ^  the 
fact  tiuit  high-voltage  power  sources  are  being  considoed  for  bodi  military  and  commocial 
space  stations,  and  the  direct  interaction  of  a  high-voltage  system  with  the  ionospheric 
plasma  could  lead  to  arcing  and  severe  power  losses. 

In  our  first  study  ^per  21),  we  conc«itrated  on  the  initial  intoaction  of  a  Inased 
sphere  with  a  magnetized,  partially-ionized  plasma  that  had  characteristics  similar  to  those 
found  at  shuttle  altitudes.  Positive  potentials  were  suddnily  applied  to  the  ^heie  and  the 
subsequoit  response  of  the  plasnui  was  modelled  via  a  numerical  solution  of  the  time- 
dqien^t,  thre^-dimensional,  nonlinear  fluid  equations  for  the  ions  and  elections  and  the 
Poisson  equation.  The  main  goal  of  the  study  was  to  determine  the  effect  of  inqract 
ionization  and  collisions  on  the  sphere-plasma  intoaction  in  the  presence  of  outgassed  or 
released  neutrals  from  the  Space  Shuttle.  Simulations  were  conducted  for  different  neutral 
species,  a  range  of  neutral  dratities,  as  well  as  different  magnetic  field  sirengdis.  The  main 
conclusions  of  the  study  are  as  follows:  (1)  When  a  hi^-voltage  sphere  is  embedded  in  a 
magnetized,  partially-ionized  plasma,  an  electron  density  torus  tends  to  form  around  the 
sphere  in  the  equatorial  plane  at  early  times.  The  torus  rotates  about  the  sphere  in  the  E  x 
B  direction.  If  the  background  neutral  gas  density  is  sufficiently  high,  the  electrons 
accelerated  by  the  strong  electric  fields  associated  witii  the  sphere  can  ionize  the  gas  by 
impact  At  a  critical  ionization  rate,  a  discharge  or  breakdown  occurs.  Collisions,  on  the 
other  hand,  act  to  decelerate  the  electrons  and  modify  the  space  charge  configuration  (i.e., 
modify  the  toroidal  shape);  (2)  For  low  neutral  densities  5  x  IQH  emr^),  collisions  and 
impact  ionization  are  not  important;  (3)  For  neutral  densities  of  the  order  of  lO^^  emr^, 
ionization  occurs  in  the  toroidal  region  where  the  electron  density  is  elevated,  but  this 
impact  ionization  does  not  appreciably  affect  the  current  collection.  In  this  case,  a  toroidal 
density  configuration  can  be  maintain^  without  a  discharge  and  die  effect  of  coUisirais  is  to 
broaden  die  torus;  (4)  For  neutral  densities  of  the  ordo:  of  emr^,  the  major  ionizatiem 
still  occurs  in  a  toroidal  region.  However,  the  ionization  causes  the  plasma  density  and 
currait  collectitMi  to  increase  npklly.  A  toroidal  discharge  is  expected  in  this  case;  (S)  For 
neutral  densities  of  the  order  of  10’^  enr^,  impact  ionization  occurs  all  around  the  sphere 
in  an  explosive  manner.  The  discharge  is  likely  to  be  spherical  in  this  case;  (6)  The  above 
results  are  for  a  magnetic  field  of  3  gauss.  If  the  magnetic  field  is  reduced  and  the  other 
parameters  are  held  fixed,  the  ionization  that  occurs  teixls  to  be  more  spherical  The  current 
collection  also  tends  to  be  higher,  and  (7)  Diffoent  molecules  have  different  ionization 
energies,  and  collision  and  ionization  cross-sections.  Of  the  four  species  studied,  Ar,  N2, 
and  O  have  very  similar  properties,  and  hence,  the  features  discussed  above  for  Ar  also 
apply  quantitatively  for  N2  and  O.  However,  barium  behaves  very  differently.  It  has  a 
larger  collision  cross-section  and  lower  ionization  energy.  Therefore,  for  the  same 
conditions,  ionization  and  discharge  processes  are  much  easier  to  trigger  in  a  Ba  plasma. 
Also,  compared  to  Ar,  O,  and  N2,  tiie  density  distribution,  as  well  as  the  discharge,  tend  to 
be  more  spherical 
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Tlie  other  two  papers  (23  &  26)  involved  high-voltage  spheres  embedded  in 
mmagnetized  plasmas.  In  one  study,  a  high  ‘positive’  voltage  was  applied  to  the  sphere 
and  the  effect  of  different  potential  rise  times  was  investigated  as  was  the  long-term 
evolution  of  the  sphere-plasma  interaction.  In  the  other  paper,  ‘negative’  voltages  were 
q)plied  to  the  sphere  and  die  subsequent  sphere-plasma  interaction  was  modelled. 

From  our  simulations  we  found  that  after  the  positive  voltage  is  applied  to  the 
sphere,  the  electrons  are  accelerated  toward  tlw  sphere  in  a  few  plasma  periods  and  an 
overshoot  oscillation  usually  follows.  The  acceleration  and  the  density  buildup  of  the 
electrons  around  the  sphere  produce  a  high  current  coUecdon  by  the  sphere.  This  high 
current  level  usually  la^  a  few  hundred  plasma  p^ods  or  a  few  ion  plasma  poiods.  Then 
the  ion  mc^on  becomes  sign«^cant  The  ions  are  accelerated  away  from  die  sphere  and  an 
ion  density  shell  forms  and  propagates  away  from  the  sphere.  The  ion  density  shell  Uocks 
die  electnm  flow  and  reduces  the  current  The  deceleration  of  die  electrons  by  die  ion  shell 
leads  to  an  electron  density  enhancen^nt  so  that  a  dcHible  layer  structure  forms.  The  double 
layer  propagates  outwards  and  eventually  slows  down  and  disiqipears.  This  process  lasts  a 
few  thousand  plasma  periods.  During  this  time,  the  current  first  decreases  to  a  minimum 
and  then  slowly  reach^  a  constant  level.  Low  frequency  oscillations  Oowo-  than  the  ion 
plasma  frequency)  nuiy  iqipear  due  to  the  overshoot  of  the  ion  density  configuration  about 
its  equilibrium  position.  However,  a  steady  state  is  eventually  reached 

In  the  steady  state,  the  ion  doisity  satisfies  die  Boltzmann  relation.  In  some  regi(»is 
near  the  sphere,  the  ions  are  completely  depleted  and  only  electrons  are  present  There  is 
also  a  region  where  the  ion  and  electron  densities  are  unequal.  The  region  with  a  net  ^ce 
charge  is  the  sheath  region.  This  region  can  be  conqiarable  to  or  larger  than  the  sphere 
radius,  depending  on  die  iqiplied  voltage.  The  sheath  thickness  estimated  from  the  electron 
and  ion  density  profiles  agrees  with  the  Langmuir-Blodgett  theory.  Far  away  from  the 
sphere,  the  plasma  is  unperturb^  Between  the  sheath  region  and  unperturbed  plasma  is 
the  presheath,  where  there  is  no  net  space  charge,  but  both  the  electron  and  ion  densities 
differ  from  the  unpoturbed  value.  The  thickness  of  the  presheath  increases  with  higher 
applied  voltages.  The  current  density  from  the  simulation  is  2  to  5  times  hi^iCT  than  that 
from  the  single  Langmuir-Blodgett  theory. 

The  method  of  applying  the  voltage  on  the  sphere  affects  die  tenqxiral  evolution.  If 
the  voltage  is  applied  according  to  an  exponential  law  with  a  rise  time  t,  only  the  election 
motion  at  very^euly  times  is  affected  if  t  ^  1.  For  t  =  0. 1 ,  the  result  is  almost  identical  to 
the  result  for  x  =  0  (a  step  function  application  of  the  voltage).  ^For  x-  1.0,  the  electron 
density  and  velocity  near  the  sphere  are  Iowa-  than  those  for  x  =  0  or  0.1  before  five 
plasma  poiods,  and  so  is  the  current  Afterwards,  the  results  are  the  same.  For  x  =  10, 
there  are  several  differences.  The  acceleration  of  the  electrons  is  much  lower.  The 
overshoot  oscillation  at  early  times  does  not  occur,  pie  ion  motion  is  affected.  The 
formation  of  the  density  shell  is  (klayed.  Fot  different  x’s,  the  system  reaches  the  steady 
state  in  slightly  different  ways.  For  a  small  x,  the  ion  density  configuration  may  pass  the 
final  equilibrium  location  and  oscillate  about  it  For  a  large  x,  such  an  overshoot  oscillation 
is  more  likely. 
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In  our  paper  dealing  with  high  ‘negative’  voltage  spheres,  simulations  were 
conducted  for  a  range  of  voltages  (-10  to  -10,000  V)  and  for  plasma  densities  from  lO^  to 
10^  cm~3.  The  temporal  evolution  of  the  plasma  was  followed  all  the  way  to  the  steady 
state.  In  all  simulations,  there  were  certain  qualitative  similarities  in  the  plasma  response. 
At  very  early  times,  the  rapid  electron  motion  away  from  the  sphere  resulted  in  an  electron 
overshoot  oscillation  (ringing).  The  frequency  of  this  oscillation  was  a  fraction  of  the 
electron  plasma  frequency  and  increased  slowly  with  the  magnitude  of  the  voltage  on  the 
sphere.  For  high  plasma  densities,  the  oscillation  was  damped  as  the  initial  ion  current 
surge  to  the  sphere  developed.  Subsequently,  the  ion  current  reached  a  peak,  decreased, 
and  then  attained  a  steady-state  level.  The  final  sheath  and  presheath  sizes  for  negative- 
voltage  spheres  were  similar  to  posidve-voltage  spheres  with  Ae  same  potential  magnitude, 
although  the  temporal  evolution  of  the  plasma  was  different  for  positive  and  negative 
spheres.  In  the  st^y  state  and  for  a  given  plasma  density,  the  ion  density  structure  in  the 
sheath  varied  with  the  applied  voltage  on  the  sphere.  As  the  voltage  increased,  an  ion 
density  hole  developed  around  the  sphere  in  the  sheath  region.  Near  the  sphere  the  ion 
density  was  relatively  high,  but  the  ion  density  decreased  with  distance,  reached  a 
minimum,  and  then  increased  to  the  background  value.  For  a  -10,(X)0  V  sphere,  the  ion 
density  near  the  sphere  was  also  four  times  larger  than  the  background  plasma  densiy . 


4.5  Magnetosphere  -Ionosphere  Coupling  via  Electric  Fields 

It  is  well  known  that  the  electric  fields,  particle  precipitation,  auroral  conductivity 
enhancements,  and  Birkeland  currents  that  couple  the  magnetosphere-ionosphere  system 
are  strongly  dependent  upon  the  direction  of  the  IMF.  When  the  IMF  is  southwaixl,  the 
Birkeland  (or  field-aligned)  currents  flow  in  the  Region  1  &  2  current  sheets,  the  F  region 
plasma  convection  exhibits  a  2-ceil  structure  with  antisunward  flow  over  the  polar  cap,  the 
conductivity  enhancements  are  confined  to  the  statistical  auroral  oval,  and  the  auroral 
electron  precipitation  is  also  confined  to  the  classical  oval.  However,  when  the  IMF  is 
northward,  the  situation  is  considerably  more  complicated  and  less  clear.  In  this  case,  an 
additional  field-aligned  current  system  occurs  in  the  polar  cap  called  the  NBZ  currents; 
plasma  convection  can  be  sunward  in  the  polar  cap  and  the  pattern  can  assume  multi-cell, 
severely  distorted  two-cell  or  turbulent  characteristics;  and  particle  precipitation  occurs  in 
the  polar  cap  that  can  be  uniform,  in  the  form  of  multiple  sun-aligned  arcs,  or  in  a  B-aurora 
configuration. 

We  have  been  and  are  continuing  to  conduct  several  theoretical  studies  in  an  effort 
to  elucidate  the  ionosphoe-magnetosphere  coupling  processes  during  northward  IMF.  One 
of  our  efforts  involved  the  use  of  our  electrodynamic  model  to  investigate  the  ‘large-scale’ 
Held-aligned  currents  that  exist  in  the  polar  cap  during  northward  IMF  (papers  24  &  34). 
With  our  electrodynamic  model,  we  solve  Ohm’s  law  and  tfie  current  continuity  equation  so 
as  to  obtain  self-consistency  between  the  ionospheric  conductivity,  field-aligned  current, 
horizontal  E-region  current,  and  convection  electric  frvld.  In  one  effort,  we  used  the  USU 
conductivity  model  and  conducted  a  systematic  study  of  the  influence  of  the  ionospheric 
conductance  tm  the  ftxm  of  the  freld-aligned  current  associated  with  the  Heppner-Maynard 
‘distorted  two-cell’  convection  pattern  (paper  24).  Our  modelling  results  indicated  that. 
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contrary  to  previous  claims,  the  NBZ  current  can  be  associated  with  the  distorted  two-ceU 
convection  pattern  for  most  of  the  ionospheric  conductivity  conditions.  We  found  that  the 
seasonal  and  auroral  activity  conditions  significantly  affea  the  iono^heric  conductivity  and 
that  the  conductivity  variations  can  influence  die  basic  features  of  the  NBZ  current  system 
associated  with  the  distorted  two-cell  convection  pattern.  Based  on  these  results,  we 
suggested  that  the  field-aligned  current  system  oteerved  by  the  MAGS  AT  satellite  might 
inqily  a  distorted  two-cell  convection  pattern,  and  that  a  four-cell  pattern  is  more  likely  to 
occur  when  the  IMF  is  due  north  or  voy  close  to  the  north. 


4.6  Magnetosphere  -  Ionosphere  Coupling  via  Ion  Outflow 

The  outflow  of  thermal  plasma  from  the  topside  ionosphov  at  high  latitudes  (i.e., 
the  polar  wind)  is  very  important  because  the  escaping  thermal  plasma  drains  the 
ionosphere  and  mass  -  loads  the  magnetosph^.  However,  it  is  difficult  to  model  the 
outfiow  because  the  polar  wind  is  collision-dominated  at  low  altitudes,  passes  through  a 
transition  region,  and  then  becomes  collisionle^  at  high  altitudes.  Over  the  years,  various 
mathemadcals  have  been  used  to  describe  the  outflow  (hydrodynamic,  collisionless  kinetic, 
generalized  transport,  etc.),  but  the  diffnent  models  have  limitations,  particularly  in  the 
transition  region  from  collision-dominated  to  collisionless  flow.  In  an  effort  to  better 
understand  the  strengths  and  limitations  of  the  various  models,  we  have  continued  our 
effort  whereby  different  mathematical  models  are  used  for  the  same  polar  wind  conditions 
so  diat  a  direct  comparison  of  results  can  be  made. 

During  the  last  year  we  completed  three  papers  dealing  with  the  validity  of  the 
various  mathematical  m^els  of  the  polar  wind  and  related  flows  (papers  20, 25,  &  32).  In 
one  study  (paper  25),  we  compared,  in  as  consistent  a  maimer  as  possible,  solutions  to  the 
bi-Maxwellian-based  16'moment  set  of  transport  equations  with  those  obtained  from  a 
semikinetic  model,  assuming  boundary  conditions  characteristic  of  both  supersonic  and 
subsonic  flows  in  the  terrestrial  polar  wind  as  well  as  supersonic  flow  in  the  solar  wind. 
For  each  case  in  which  transport  and  senukinetic  solutions  were  compared,  three  separate 
semikinetic  solutions  were  generated.  These  three  semikinetic  solutions  assumed  the 
particle  distribution  function  at  the  baropause  to  be  an  isotropic  Maxwellian,  a  bi- 
Maxwellian,  and  a  bi-Maxwellian-based  16-moment  expansion  with  zero  stress, 
respectively.  Our  study  demonstrated  several  important  points;  (1)  For  supersonic 
“collisionless”  flows,  the  16-moment  transport  theory  and  the  semikinetic  model  assuming 
a  16-moment  distribution  at  the  baropause  are  almost  identical  in  their  predictions,  even  for 
the  higher-order  moments  (parallel  and  perpendicular  heat  flows).  This  is  true  for  both 
polar  and  solar  wind  conditions.  (2)  The  semikinetic  solutions  assuming  either  a 
Maxwellian  or  a  pure  bi-Maxwellian  at  the  baropause  also  show  octremely  close  agreement 
with  the  transport  results  for  the  lower-order  moments  (density,  drift  velocity,  and  parallel 
and  perpendicular  temperatures),  but  are  less  accurate  in  their  heat  flow  predictions.  (3) 
The  nearly  precise  agreement  between  the  16-moment  transport  solutions  and  the 
semikinetic  solutions  with  a  16-moment  distribution  at  the  boundary,  which  implicitly 
contain  the  full  hierarchy  of  moment  equations,  indicates  that  moments  higher  than  heat 
flow  (flow  of  parallel  and  perpendicular  thermal  energy)  are  not  needed  to  describe  the 
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steady-state  polar  and  solar  wind  cases  considered  in  this  study.  (4)  Because  of  its 
underlying  assumptions,  the  semikinetic  model  is  unable  to  properly  describe  subsonic  H*** 
flows.  Therefore,  a  comparison  of  semikinetic  and  transport  models  for  subsot.,  j  flow 
conditions  must  await  future  advances  in  the  kinetic  theory.  (S)  Collisions  are  clearly  of 
importance  in  determining  the  thermal  and  heat  flow  structure  of  the  solar  wind.  The  16- 
moment  transport  model,  which  incorporates  the  effects  of  Coulomb  collisions,  yields 
temperature  anisotropies  at  1  AU  that  are  in  agreement  with  measurements,  while  the 
serrdkinetic  model,  Wch  is  collisionless,  does  not 
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5.  FUTURE  DIRECTION 

During  the  third  year  of  the  program,  we  will  continue  our  collaborative  projects 
widi  the  various  scientists  at  the  I%illips  Laboratory  in  Bedford,  Massachusetts.  We  will 
also  continue  research  in  die  following  general  areas: 

1 .  Ionos[dieie  Stmcture  Modelling 

2.  Operational  lonosphoic  Models 

3.  Model^DataGcMiqmistxis 

4.  OTH  Radar  Studies 

5.  Direct  Magnetosidieric  Monitoring  hidex 

6.  Plasma  Convection  Studies 

7.  Magnetosphere-I(Hio^ere  Coqiling 

8.  Thennospheric  General  Grculadmi  Modelling 

9.  Plasma^here  Dynamics 

1 0.  Spacecraft-Environment  InteracticHi  Studies 
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1 .  American  Geophysical  Union  Fall  Meeting 
San  Francisco,  Qdifomia 

12/4- 12/8/89 

Schunk,  Sojka,  Yang,  Demars,  Ma,  Rasmussen,  and  Thionann  attended  meeting  and 
13  papers  were  presented. 

2.  Space  Model  Devek^Mimit  Review  Meeting 
Fetoson  Air  Force  Base,  Colorado  Springs,  Colorado 
1/29/90 

Schunk  and  Sojka  attended  meeting  and  a  joint  paper  was  presented. 

3.  APGLMeedng 
Boston,  Massachusetts 
3/12-3/13/90 

Schunk  and  Sojka  attended  meeting  to  discuss  collaborative  projects  with  D. 
Anderson,  F.  Marcos  and  N.  Maynard. 

4.  SUNDIAL  Meeting 
New  Oiieans,  Lot^iana 
4/9-4/12/90 

Schunk  attended  meeting  and  presented  an  invited  pq>er. 

5.  I(»oiq)heric  Specification  Model  Quarteriy  Review 
AFGL;  Bosttm,  Massachusetts 

5/14/90 

Schunk  and  Sojka  presented  talks. 

6.  CEDAR  Woricshop 
Boulder,  Colorado 
6/15/90 

Sojka  was  a  panel  member  involving  High  Latitude  Plasma  Structures. 

7.  Solar-Terrestrial  Physics  Symposium 
The  Hague,  The  Neiheilands 

6/25  -  6/30/90 

Schunk  presented  an  invited  review  talk. 

8 .  COSPAR  Intnnational  Meeting 
The  Hague,  The  Netherlands 
7/2-7/6/90 

Schunk,  Sojka  and  Thiemann  attraded  and  8  piq)ers  wnc  presented. 

9.  Camtmdge  Worktop  on  Tiansptnt  in  Geoplasmas 
CamMdge,  Massachusetts 

7/16-7/20/90 

Barakat  attended  and  presented  2  papers. 

10.  Gordon  Research  Conference 
Plymouth,  New  Hampshire 
7/30-8/3/90 

Schunk  {vesoited  an  invited  talk. 
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11.  AFOSR  Siqyported  Research 
Boulder,  Colorado 
8/30-9/1/90 

D.  Payne  visited  NCAR  to  discuss  the  development  of  an  auroral  deposition  code 
with  Dirk  Lummctzheiia 

12.  AFGL  Meeting 
Boston,  Massachusetts 
9/3-9/WW 

Schunk  and  Sojka  discussed  collaborative  efforts  with  H.  Carlson,  W.  Doiig,  R. 
Daniell,  D.  Hardy,  S.  Gussenhoven,  and  J.  Whalen. 

13.  Ioooq)hm  Spedficadon  Model  Quarterly  Review 
Scott  Air  Force  Base,  Illinois 
9/17-9/18/90 

Sojka  attended  meeting  and  presented  talk. 

14.  Piene-LouisBldly  visit  to  USU 
Logan,  Utah 

O^ober  1990 

To  work  on  a  project  widi  R.W.  Schunk. 

15.  AFOSR  Supported  Research 
Btmider,  Colorado 

11/13  -  1V16/90 

D.  Payne  visited  NCAR  to  compare  auroral  deposition  codes  with  Dirk 
Lumnerdieita 

16.  H.  Thiemann  visit  to  USU 
Logan,  Utah 

November-  Deoentber  1990 

To  work  on  several  projects  with  RW.  Schunk. 

17.  American  Geophyi^  Union  Fall  Meeting 
San  Francisco,  California 
12/2-12/7/90 

Demars,  lOioyloo,  Ma,  Thiemann,  and  Yang  attended  meeting  and  presented  piq>ers. 

18.  AFOSR  Sunrarted  Research 
Hanscom  Air  Force  Base,  Massachusetts 
12/11-12/19/90 

P.  Qtrone  visited  die  Get^ysics  Laboratory  to  cdlect  data  and  discuss  research. 

19.  Conference  on  C^  Systems 
Monterey,  Califtmiia 
2/9-2/17/91 

Schunk  attended  meeting  and  presented  two  papers. 

20.  lonoqdieric  Specification  Model  Quarterly  Review 
Coloiw  Springs,  Colorado 

3/4-3/7/91 

J  J.  Sojka  attended  meeting. 
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21.  PRIMOWoikshop 
Bottldo:,  Colorado 
6/20-6/21/91 

R.W.  Schunk  and  J.J.  Sojka  attended  this  workshop  which  was  organized  by  D. 
Anderson  of  Phillips  Laboratoiy. 

22.  AFOSR  Supported  Research 
Boston,  Massachusetts 
7/18  -  7/20/91 

J.J.  Sojka  visited  the  Geophysics  Laboratory  to  discuss  collaborative  research 
projects. 

23.  lUGG  General  Scientific  AssemUy 
Vienna,  Austria 
8/16-8/22/91 

H.G.  Detnars  attended  meeting  and  presented  paper. 

24.  AFOSR  Siqypoited  Research 
Fairbanks,  Alaska 

9/17  -  9/19/91 

L»  Zhu  visited  Joe  Kan  at  the  University  of  Alaska  to  discuss  collaborative  research 
projects. 


